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ABSTRACT 

GO.253+0.016 is a molecular clump that appears to be on the verge of forming 
a high mass, Arches-like cluster. Here we present new ALMA observations of its 
small-scale (M).07 pc) 3 mm dust continuum and molecular line emission. The 
data reveal a complex network of emission features, the morphology of which 
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ranges from small, compact regions to extended, filamentary structures that are 
seen in both emission and absorption. The dust column density is well traced by 
molecules with higher excitation energies and critical densities, consistent with 
a clump that has a denser interior. A statistical analysis supports the idea that 
turbulence shapes the observed gas structure within GO.253+0.016. We find 
a clear break in the turbulent power spectrum derived from the optically thin 
dust continuum emission at a spatial scale of WJ.l pc, which may correspond 
to the spatial scale at which gravity has overcome the thermal pressure. We 
suggest that GO.253+0.016 is on the verge of forming a cluster from hierarchical, 
filamentary structures that arise from a highly turbulent medium. Although the 
stellar distribution within Arches-like clusters is compact, centrally condensed 
and smooth, the observed gas distribution within GO.253+0.016 is extended, with 
no high-mass central concentration, and has a complex, hierarchical structure. 

If this clump gives rise to a high-mass cluster and its stars are formed from this 
initially hierarchical gas structure, then the resulting cluster must evolve into a 
centrally condensed structure via a dynamical process. 

Subject headings: dust, extinction-stars:formation-ISM:clouds-infrared:ISM-radio 
lines: ISM 


1. Introduction 


Young Massive Clusters (YMCs, e.g., Arches, Quintuplet, Westerlund 1, RSCG, GLIMPSE- 
COl: Figer et al. 1999; Clark et al.||2005 ; Figer et al. 2006; Davies et al. 2011) are gravita¬ 


tionally bound stellar systems with masses >10 4 M Q and ages < 100 Myr (Portegies Zwart 


et al. 2010). Despite their small numbers and short lifetimes, they appear to be the ‘missing 
link’ between open clusters and globular clusters (e.g., Elmegreen fc Efremov||1997 ; Bastian 


2008; Kruijssen 2014). As such, understanding their formation may provide the link to un¬ 


derstanding globular cluster formation. If YMCs and even globular clusters form in a similar 
way to open clusters, then our knowledge of nearby stellar systems may potentially be used 
as a framework to understanding star formation across cosmic time, back to globular cluster 
formation ~ 10 Gyr ago. 

An understanding of the formation of clusters requires observations of their natal dust 
and gas well before the onset of star formation. In recent surveys of dense cluster-forming 
molecular clumps, one object, GO.253+0.016, stands out as extreme (Longmore et al. 2012). 
Located at a distance of ~8.5kpc, GO.253+0.016 lies ~100pc from the Galactic centre 
(Molinari et al. 2011). With a dust temperature of ~20 K, effective radius of ~2.9pc, 
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average volume density of >10 4 cm 3 , mass of ~10 5 M 0 , (Lis fc Carlstrom|l994; Lis & Menten 


1998; Lis et al. 2001) and devoid of prevalent star-formation, GO.253+0.016 has exactly the 
properties expected for the precursor to a high-mass, Arches-like cluster. Despite being 
identified over a decade ago, its significance as a precursor to a high-mass cluster has only 
recently been investigated in detail (Longmore et al.|2012; Kauffmann et al.|2013; Rathborne 


et al. 2014a; Johnston et al. 2014). 


Recent work shows that the volume density threshold predicted by theoretical models 
of turbulence, which depend on the gas mean volume density and 1-dimensional turbulent 
Mach number, is ~10 4 cm -3 for the solar neighbhourhood and is ~10 8 cm -3 for the high- 
pressure inner 200 pc environment of the Central Molecular Zone (CMZ; Kruijssen et al. 


2014b; Rathborne et al. 2014b). The lack of current star formation within GO.253+0.016 is 


consistent with its location in the CMZ and this environmentally dependent volume density 
threshold for star formation, which is orders of magnitude higher than that derived for the 


solar neighbourhood (Rathborne et al. 2014b). 


Because GO.253+0.016 is close to virial equilibrium, has sufficient mass to form a YMC 
through direct collapse, and shows structure on ~ 0.3pc scales (Longmore et al. 2012), it 


provides an ideal template to investigate whether or not YMCs form through hierarchical 
fragmentation, in a scaled up version of open cluster formation. Detailed studies of clumps 
like GO.253+0.016 have the potential to reveal the gas structure in protoclusters before star 
formation has turned on and significantly disrupted its natal material. If YMCs like the 
Arches are the missing link between globular clusters and more typical open clusters, the 
detailed study of protoclusters like GO.253+0.016 may reveal the initial conditions for star 
formation not only within high-mass clusters, but potentially for all clusters. 


Both theoretical studies (Elmegreen 2008; Kruijssen 2012) and simulations (Bonnell 


et al. 2008; Krumholz et al. 2012; Dale et al. 2013) suggest that high-mass clusters form 


through the hierarchical merging of overdensities in the interstellar medium. While some of 


this hierarchical structure may also be evident in the initial stellar distribution (Maschberger 


et al. 2010; Hopkins 2012; Parker et al. |2014| ), it can be easily erased over long time-scales, 
either via the structure’s dispersal (if unbound) or relaxation (if bound). In either case, 
these processes are expected to give rise to spherically symmetric stellar clusters and a 
background field star population (e.g. Lada & Lada 2003). Notwithstanding this theoretical 
picture, numerical studies of the early phase of cluster evolution often model the gas as a 


spherically-symmetric background potential (e.g. Goodwin & Bastian 2006; Baumgardt &; 
Kroupa||2007 ), although notable exceptions exist (e.g. Moeckel fe Bate||2010 ; Kruijssen et al. 
2012; Parker & Dale 2013). This disparity between cluster formation theory, which suggests 
a highly structured initial state, and the early evolutionary modelling, which usually posit 
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a smooth spherical initial state can be resolved by characterising the observed gas structure 
within high-mass protoclusters before star formation has begun. 

The outcome of the high-mass cluster formation process is well-known. Observations 
show that YMCs are centrally condensed with > 10 4 M 0 of stars in their central ~ 1 pc. 
Thus, the puzzle is then: is this central condensation of stars set by the initial gas structure 
or do they relax into this small, centrally concentrated region over time? 

To answer this question and to distinguish between the models, requires observations 
of protocluster clumps with a significant improvement in angular resolution, dynamic range, 
and sensitivity compared to those obtained to date. Fortunately, the new Atacama Large 
Millimeter/Sub-millimeter Array (ALMA) provides the necessary improvements and allows 
us, for the first time, to image protoclusters across the Galaxy with comparable physical 
resolution (i.e., < 0.1 pc) to clouds in the nearby solar neighbourhood. 

Our goal is to determine the location, mass, and kinematics of the small-scale fragments 
within GO.253+0.016. The new ALMA data reveal the internal structure and kinematics in 
this high-mass protocluster down to spatial scales of 0.07 pc (Fig. [Tj) . As such, these data 
can be used to determine how the material is distributed within this protocluster. 


2. ALMA observations 

We imaged the 3 mm (90 GHz) dust continuum and molecular line emission across 
GO.253+0.016 using 25 antennas as part of ALMA’s Early Science Cycle 0. Because ALMA’s 
field of view at this wavelength is ~ 1.2’ and the cloud is large on the sky (3’ x 1’), a 13- 
pointing mosaic was obtained to cover its full extent. 

Using ALMA’s Band 3, the correlator was configured to use 4 spectral windows in dual 
polarization mode centred at 87.2, 89.1, 99.1 and 101.1 GHz each with 1875 MHz bandwidth 
and 488 kHz (1.4—1.7 km s _1 ) channel spacing. With this setup we obtained the 3 mm 
dust continuum emission and also covered a number of excellent tracers of dense molecular 
and shocked gas, in particular: HCN, HCO+, C 2 H, SiO, SO, H 13 CO+, H 13 CN, HN 13 C , 
HCC 13 CN, HNCO, H 2 CS CH 3 SH, NH 2 CHO, CH 3 CHO, NH 2 CN, and HC 5 N (for brevity in 
the text, we refer to the molecular transitions simply by the name of the molecule, e.g., HCN 
v=0 J=l—0 will be called ‘HCN’; see Table [l] for the quantum numbers of the transitions, 
frequency, critical densities, and excitation energies). 

To build up the signal to noise and improve the uu-covcragc, GO.253+0.016 was imaged 
on six different occasions over the period 29 July - 1 August 2012. We used the ALMA early 
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science ‘extended configuration’; the projected baselines range from 13-455 m. All data 
reduction was performed using the CASA and Miriad software packages. Each individual 
dataset was independently calibrated (by our ALMA support scientist) before being merged 
together. We performed the cleaning and imaging to produce the final data products. Be¬ 
cause the emission from adjacent 1.7km s -1 spectral channels are correlated, the final data 
cubes were resampled to have a velocity width of 3.4 km s -1 per channel. The ALMA data 
have an angular resolution of ~1.7" (~0.07 pc), with a la continuum sensitivity of ~ 25//.Jy 
beam -1 and line brightness sensitivity of 1.0 mJy beam -1 per 3.4 km s -1 channel. 


3. Combination of the ALMA interferometric data with single dish data 


Because the interferometric ALMA data are insensitive to emission on the largest angu¬ 
lar scales (> 1.2'), where possible we used single-dish data to fill in the missing zero-spacing 
information. Given that the emission from the clump is extended (3'xl') and the emission 
from the dust continuum and many of the molecular transitions on spatial scales > 1.2' 
significant, it was critical that the zero-spacing information was included to properly sample 
the wu-plane (e.g., 


Stanimirovic 2002 


The combination of the ALMA data with the zero-spacing information was performed 
in the Fourier domain and was achieved using the corresponding single-dish data as a model 
during the CLEANing process. The inclusion of the single-dish data in this manner is 
superior to using alternatives task such as ‘feathering’ (see Ott 2014 for details). For this 
technique to be applicable there must be an overlapping region of spatial frequencies between 
the two datasets. Fortunately, for both the continuum and line datasets, existing single-dish 
data provided the necessary coverage for inclusion with the ALMA data. 


Ideally, the single dish data should have the same signal-to-noise as the interferometric 
data that they are being combined with. Although the signal-to-noise of the single dish data 
is lower than that of the ALMA data, it is high enough that the gain from recovering the 
zero-spacing information outweighs the increased noise that they may introduce. 


To recover the zero-spacing information for the ALMA line cubes (HCN, HCO + , C 2 H, 
SiO, HNCO), we used data from the MALT90 survey which were obtained using the Mopra 
22m telescope (Foster et al. 2011; Jackson et al. 2013; Rathborne et al. 2014a). To convert 
these data from T^ to the flux F v in Jy beam -1 , we used a scaling factor of 23 K/Jybeam -1 , 
calculated via the equation T^= 1.36 A 2 F v /( r]0 2 ), using a FWHM beam size, 9, of 38", 
wavelength, A, of 3mm, and a main beam efficiency, 77, of 0.49 (Ladd et al. 2005). 


For the continuum image, since we do not have a direct measurement of the emission at 
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3 mm on the larger angular scales, it was necessary to estimate the flux at 3 mm. To achieve 
this we scaled an image of the dust continuum emission at 500 /mi ( Herschel , 33" angular 
resolution) to what is expected at 3 mm assuming that the flux scales like u (2+B K where /3 
was assumed to be 1.2. Herschel data were used for this purpose as these data provide the 
most reliable recovery of the large scale emission. We choose to scale the 500 /mi Herschel 
image as it is the closest in wavelength to 3 mm, minimising the effect of the assumption of 
j3 (see Rathborne et al. 2014b for more details). Figure [2] compares the single dish, ALMA- 
only, and combined (ALMA and single dish) images for both the continuum and the HNCO 
integrated intensity. For the 3 mm dust continuum emission, the ALMA-only image recovers 
~18% of the total flux. The images clearly show the significant improvement in angular 
resolution provided by ALMA and the recovery of the emission on the large spatial scales 
provided by the inclusion of the zero-spacing information. 


4. The small-scale structure within a cluster forming clump revealed 

4.1. Continuum emission 


Figure [I] shows a two colour image of GO. 253+0.016 combining infrared images from the 
Spitzer Space Telescope (24 /im, ~6" angular resolution, shown in cyan) with the new mil¬ 
limetre image from ALMA (3mm, ~1.7" angular resolution, shown in red). Contours show 
single dish dust continuum emission (JCMT 450/im, ~ 7.5" angular resolution). Identified 
as an infrared dark cloud (IRDC), GO.253+0.016 is clearly seen as a prominent extinc¬ 
tion feature in mid-IR images. Because it is so cold, its thermal dust emission peaks at 
millimeter/sub-millimeter wavelengths. Moreover, since the 3 mm dust continuum emission 
is optically thin, internal structure is easily revealed by the ALMA images. While the over¬ 
all morphology of the dust continuum emission revealed by ALMA matches well the mid-IR 
extinction and the large-scale dust continuum emission provided by the single dish image, 
the ALMA image shows significant structure on the small-scale. The emission is neither 
smooth nor evenly distributed across the clump; it is highly fragmented and there is a clear 
concentration of emission toward its centre. These small-scale fragments are presumably the 
cores from which stars may form. 


Because the 3 mm continuum emission is optically thin , given standard assumptions 
about the dust emissivity, it can be used to determine the column density and, hence, mass 
of its small-scale fragments. Recent observations from Herschel of the 170-500 /un dust 
continuum emission, show that the dust temperature (T*^) within GO.253+0.016 decreases 
smoothly from 27 K on its outer edges to 19 K in its interior (Longmore et al. 2012). External 
to GO.253+0.016, the dust temperatures are considerably higher, >35 K. Within the cloud 








and on the spatial scales probed by the Herschel data (~33", ~1 pc), no small volumes of 
heated dust are apparent. All of the millimetre emission detected by ALMA falls within the 
region for which was measured to be < 22 K. Thus, we assume that the temperature 

of the dust detected by ALMA ranges from 19-22 K and use a value of 20 ± 1K for all 
calculations. 


Assuming this dust temperature, a standard value of the gas to dust mass ratio of 100, a 
dust absorption coefficient (n 3 mm ) of 0.27 cm 2 g -1 (using Ki. 2mm = 0.8 cm 2 g -1 , and k oc 
Chen et al.|2008 ; Qssenkopf fe Henning|1994 ,) and a f3 of 1.2 ± 0.1 , the intensity of the 3 mm 


emission can be converted into the H 2 column density, N(H 2 ), by multiplying the measured 
fluxes by a single conversion factor of 1.9 x 10 23 cm -2 /mJy. We assume that T dust is constant 
across the cloud. For the small range of measured dust temperatures in GO.253+0.016, this 
approximation is well justified. Given these assumptions, the uncertainties for T dust and d 
introduce an uncertainty of 10% for the column density, volume density, mass, and virial 
ratio. 


The column densities across GO.253+0.016 measured from the ALMA-only image range 
from 0.5 - 82 x 10 22 cm -2 , with a mean value of 1.0 x 10 22 cm -2 . Since this image detects only 
the highest contrast peaks within the cloud and does not include the large-scale emission, 
we also determine the column density distribution using the combined (ALMA+SD) image 
which is more representative of the total column of material associated with GO.253+0.016. 
When including the emission on all spatial scales, we find that the mean column density is ~ 
cm -2 and an order of magnitude higher compared to molecular clouds in the solar 


8.6 x 10 22 —“ 2 


neighbourhood (see Rathborne et al. 2014b for a detailed discussion of the column density 
distribution across GO.253+0.016 and its comparison to solar neighbourhood clouds). 


4.2. Molecular line emission: spectra, moment maps, channel maps, position 

velocity diagrams 

The 3 mm (90 GHz) spectrum is rich in molecular lines that span a range in their exci¬ 
tation energies and critical densities (see Table [l|. As such, these observations also provide 
data cubes from many molecular species and reveal a wealth of information about the gas 
kinematics and chemistry within GO.253+0.016. Figure [3] shows, along with the dust contin¬ 
uum image, the integrated intensity images for each of the detected molecular transitions. 
Clearly seen are the differences in the morphology of the emission from the various tracers, 
indicating that the chemistry and excitation conditions vary considerably within the clump. 

To quantify the differences in the morphology of the integrated intensity for each molec- 
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ular transition compared to the others and the dust continuum emission, we have determined 
the 2-D cross-correlation function for all combinations of image pairs. Figure [4] shows the 
derived correlation coefficients for each transition with respect to all others (the lower half 
of the figure shows the correlation coefficients via a colour scale, the upper half of the figure 
shows the derived values). For each particular tracer (labelled along the diagonal) the molec¬ 
ular transition for which it has the highest correlation coefficient can be found by searching 
for the darkest square either in the row to its left or in the column below (the corresponding 
values are mirrored in the upper portion of the figure). 

We find that the morphology of the dust continuum emission is not highly correlated 
with any gas tracer, but is most correlated with the morphology of the emission from the 
molecules with high excitation energies (NH 2 CHO in particular). In addition, we find that 
the morphology of the emission from HCN and HCO + are highly correlated with each other 
and little else, that the shock tracers (SiO and SO) are correlated with each other, and that 
the dense gas tracers with high excitation energies are well correlated with each other. 

Figure [5] shows example spectra from the detected molecular transitions at several lo¬ 
cations across the clump (spectra from a single pixel were extracted at the positions marked 
on Fig. [3] and listed in Table [ 2 ]). The spectra clearly show the complexity of the emission 
and the differences in intensity and velocity detected from the various tracers. Molecular 
species with lower excitation temperatures (E u /k ~4K, e.g., HCN, HCO + , C 2 H) show emis¬ 
sion over a wider range in velocity (~ 10-15km s -1 ) compared to the isotolopogues (e.g., 
H 13 CN, H 13 CO + , HN 13 C , and HCC 13 CN; ~3 -10 km s _1 ) ) and those with higher excitation 
temperatures (E Jk >10 K, e.g., CH 3 SH, NH 2 CHO, CH 3 CHO, H 2 CS, NH 2 CN, and HC 5 N; 
~ 3 km s _1 ). The shocked gas (traced by SiO and SO) also shows emission over a wide range 
in velocity towards most locations in the clump (~ 3-20 km s -1 ) . 


For each molecular transition we characterise the emission via a moment analysis using 

All moments were calculated over the velocity range 


(1995). 


the definitions of Sault et al 
— 10 km s" 1 <V L5i? <70km s _1 (this range is delineated by dotted vertical lines on the 
spectra). The moment analysis allows us to characterise the emission at every point across 
the map and to easily compare the emission from the various molecular transitions to each 
other and to the dust continuum emission. 


We compute the zeroth moment (M 0 , or integrated intensity), the first moment (Mi, 
the intensity weighted velocity field), and the second moment (M 2 , the intensity weighted 
velocity dispersion, a). These are defined as 


M 0 = TX(v)dv, Mi = 


f t a( v )v dv 

jT* A (v)dv ' 


M 2 = 


1 f T* a (v)(v - M x ) 2 dv 








10 


where T^(u) is the measured brightness at a given velocity v. The peak intensity of the emis¬ 
sion was also determined via a three-point quadratic fit to each spectrum (defined as moment 
—2; Sault et al. 1995 ). Note that C 2 H has well separated hyperfine components which will 
affect the first and second moment determinations (HCN also has hyperfine components, but 
their separation is less than the observed linewidth for HCN toward GO.253+0.016). 


To include only the highest signal-to-noise emission in the moment analysis, the peak 
intensity, intensity weighted velocity field, and the intensity weighted velocity dispersion 
were only computed at positions across the clump where the emission was above an intensity 
threshold. For the emission from the molecular species that are the brightest and include 
the large-scale emission from the single dish data (i.e., HCN, HCO + , C 2 H, SiO, and HNCO) 
we used a threshold in the integrated intensity of > 20% of the peak. Since the noise level is 
different and the extended emission is included in these cases, this threshold was chosen so 
that the moments were calculated only when the emission was significantly detected above 
the noise. The choice of this threshold was visually checked (the region over which the 
moments were calculated are clearly outlined by the edges of the moment maps) to ensure 
that it was appropriate for the moment calculation. For the fainter lines with a lower signal 
to noise that do not include the large-scale emission, we used a threshold in the peak intensity 
of > 3.8 mJybeam 1 ( 5a). 


At locations within the clump where the emission arises from well separated velocity fea¬ 
tures within the range —10 km s -1 <V lsr < 70 km s -1 , the first moment (intensity weighted 
velocity field) will result in an intensity weighted average of their individual velocities. In 
addition, the second moment (intensity weighted velocity dispersion) will clearly be larger 
than the velocity dispersion of the individual velocity components. Thus, while this analysis 
may not describe well every individual velocity component, it does provide a method by 
which to easily characterise and compare the emission from the various tracers and the dust 
continuum across the clump. 


Figure [6] shows the integrated intensity image (zeroth moment, M 0 ), peak intensity im¬ 
age, intensity weighted velocity field (first moment, Mi), and also the intensity weighted ve¬ 
locity dispersion (second moment, M 2 ) derived from the HCN emission toward GO.253+0.016 
(all other lines are included on-line in Figs. 15-30). Included on the latter three images are 
contours of the peak intensity: these highlight the small regions within the clump that are 
brightest and are included to aid in the comparison to the velocity field and dispersion im¬ 
ages. All images showing the velocity field and dispersions use the same range of values for 
the colour coding (i.e., from —20 to 60km s -1 and from 0 to 20km s -1 respectively). This 
common colour coding for the different molecular transitions allows easy comparison and 
readily highlights differences between them. 










11 


The images of the velocity field clearly show the well known large-scale velocity gradient 
across the clump; most tracers of the dense gas show blue-shifted emission towards more 
positive declinations and red-shifted emission towards more negative declinations (i.e., SiO, 
SO, HNCO, and H 2 CS). Some of the more complex molecules (NH 2 CHO, CH 3 CHO, and 
H 2 CS) are also consistent with this trend. In contrast to the molecular transitions with lower 
excitation energies, the transitions from the more complex molecules also reveal emission 
toward the clump’s centre at ~20km s -1 (coloured green in these images). 


The obvious exception to these trends is seen in the emission from both HCN and 
HCO + , which are clearly red-shifted across the whole clump and with respect to all other 
tracers. The trend that these optically thick gas tracers are systematically red-shifted with 
respect to the optically thin and hot gas tracers, indicates that GO.253+0.016 exhibits radial 
motions (either expanding or contracting, depending on the relative excitation temperatures 
of the inner and outer portions of the cloud, see Rathborne et al. 2014a for an in-depth 
discussion). The trend of red-shifted HCN and HCO + emission seen in the single-dish data 
remains in the interferometric data even when the single dish data are not included. Thus, 
radial motions are evident on all spatial scales. 


Figure [T] shows the channel maps for HCN (all other channel maps are included on-line 
in Figs. [3l] - [46j ). The panels show the emission from each spectral channel (the correspond¬ 
ing velocity is labelled in the top left corner of each panel) toward positions within the 
image where the emission was above the previously-defined intensity threshold. The channel 
maps also show well the clump’s large-scale velocity gradient. On the smaller scales, they 
reveal a complex network of emission features with a complicated velocity structure: there 
is emission on all spatial scales, the morphology of which ranges from small, compact re¬ 
gions to extended, filamentary structures. Toward many of the filamentary features we find 
homogeneous velocity fields, indicating that they are in fact physically coherent structures. 


Surprisingly, we find filamentary features in both emission and in absorption. The 
absorption filaments are evident in the HCN and HCO + integrated intensity images and 
channel maps (i.e., Figs. [6, 15 and Figs. [71 [M] respectively). The most intriguing absorption 
filaments are those seen in HCO + that extend for more than 20 km s -1 in velocity. A detailed 
analysis of these features is discussed in Bally et al. (2014). 


Position-velocity diagrams were also generated for each of the detected molecular tran¬ 
sitions (Fig. [8]). The major axis of the clump is defined along the 0-offset position in Right 
Ascension. The position-velocity diagrams were made by averaging the emission over the 
clump’s minor axis (i.e., in Right Ascension) and show the emission along the major axis 
(Declination) as a function of velocity. These diagrams also show the complexity in the 
emission and the systematic red-shift of the HCO + and HCN emission. 
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5. Statistical characterisation of the emission 

Given the complexity of the clump’s structure and velocity field, we have employed 
several statistical methods in order to characterise the observed emission from the dust 
continuum and the various molecular transitions. Our goal is to answer the questions: How 
correlated is the emission from the various tracers with respect to one another? Which 
molecular line, if any, traces well the dust column density? Does turbulence set the observed 
structure? What is the structure of the gas; is it filamentary and hierarchically distributed 
or smooth and centrally condensed? Answers to these questions will reveal both the physical 
conditions within GO.253+0.016 and the mechanisms that are shaping its internal structure. 


5.1. Molecular ‘families’ tell the same story 

The combination of molecular transitions covered by these observations trace distinct 
chemical and physical conditions. For comparison, we group the molecules into four ‘families’: 
(1) those that have ‘low’ excitation energies E u /A; ~ 4K (i.e., HCN, HCO + , C 2 H), (2) 
those that are tracers of shocked gas (i.e., SiO and SO), (3) the isotopologues (i.e., H 13 CN, 
H 13 CO + , HN 13 C , and HCC 13 CN), and (4) those that have ‘high’ excitation energies E u /k 
> 10K (i.e., HNCO, CH 3 SH, NH 2 CHO, CH 3 CHO, H 2 CS, NH 2 CN, and HC 5 N). 

Figures 10 compare the derived parameters from the moment analysis to the dust 
column density. These plots were generated by extracting the values from the moment maps 
and column density image at every pixel where there were non-zero values in both images . 
Because the number of plotted points within each panel is high, the data were binned so that 
a contour map of the distribution could be made. Contours show the overall distribution and 
location of the most common values (the contour levels are shown at 10 to 90% in steps of 10% 
of the peak number; the dotted contour outlines the 1% level). From top to bottom the rows 
show the normalised integrated intensities, normalised peak intensities, intensity weighted 
velocity field, and intensity weighted velocity dispersion as a function of dust column density. 
For ease of comparison, the molecular transitions are grouped according to the molecular 
‘families’. Within each group the molecular transitions are ordered (from left to right) by 
increasing excitation energies (see Table [l] for the specific values). 

Figure [9] shows the plots for HCN and HCO + . For comparison, we show the same 
plots made including and excluding the large-scale emission (ALMA+SD and ALMA only; 
left and right respectively) to show the effect of its inclusion when calculating the intensity 
weighted moments. As expected, when including the large-scale emission, we find a clear 
correspondence between HCN and HCO + and the lower dust column density material (i.e., 
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at column densities < 5 x 10 22 cm -2 ). These molecules clearly show an absence of blue- 
shifted emission, compared to most others (compare to Figs. [lO] and [IT]). When the large-scale 
emission is included, the measured intensity weighted velocity dispersion is more centrally 
peaked (because the emission is widespread, with a broad line profile, the ‘intensity weighting’ 
will tend to result in a similar calculated value for the intensity weighted velocity despite 
the emission being detected over a broad range). Moreover, a clear trend is evident: the 
measured velocity dispersions decrease as the dust column density increases. This trend 
is not evident in the plots made with the ALMA-only data because, in this case, there is 
no emission from the molecular transitions at the same location as the low column density 
emission, which is required to see the trend. 

The shock tracers (SiO and SO) show bright and widespread emission in both position 
and velocity across GO.253+0.016 (Fig. [lO]). Moreover, their measured velocity dispersions 
extend to ~ 20 km s -1 , which is higher than most other tracers and expected since SiO and 
SO trace gas that is stirred up in regions of shocks. When including the large-scale emission, 
SiO shows the opposite trend in its velocity dispersion compared to HCN and HCO + : its 
velocity dispersion increases as the dust column density increases, consistent with increased 
shocks in the higher density material. For SO, while its typical measured dispersion is low (~ 
2 km s -1 ), there is a significant tail in this distribution up to dispersions of ~ 20 km s -1 . In 
contrast to SiO, there is no clear trend for an increasing SO velocity dispersion with increased 
dust column density: SO shows low velocity dispersions regardless of the column density. 
The differences in the velocity dispersions measured for SiO and SO could be understood 
in terms of the differences in their excitation energies and critical densities : while SO has 
a higher excitation energy compared to SiO, it has a critical density which is an order of 
magnitude lower. Thus these species may in fact be tracing different material within the 
clump. 


There are many filamentary features found within GO.253+0.016 that are detected in 
SiO and SO, that have high velocity dispersions, and large steps in velocity over small spatial 


regions (i.e., linear features seen in the moment maps in Figs. 18 and 19). Several of these 
are also evident in the channel maps of H 13 CN (at LSR velocities of ~ 10-17 km s -1 ; see 
Figs. [34j [35j |36j) . They likely correspond to shock fronts within the clump: their details will 
be presented in a future paper. 

The isotopologues and the molecular transitions with higher excitation energies (e.g., H 13 CN, 
H 13 CO+, HN 13 C , HCC 13 CN, HNCO, CH 3 SH, NH 2 CHO, CH 3 CHO, H 2 CS, NH 2 CN, and 


HC 5 N; Figs. 10 and 11), are detected over a much smaller range in dust column densities 


and are more sharply peaked at, or slightly above, the mean value of the dust column density 
(~ 8.6 xlO 22 cm -3 , shown as the vertical dotted line) compared to the molecules with lower 
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excitation energie^j] Their integrated intensity distribution appears bi-modal, indicating the 
presence of two distinct integrated intensity values. Their measured velocities span from 0 
- 50 km s -1 ; this velocity distribution is broader compared to that derived from HCN and 

HCO+0 


While the peak intensities for the isotopologues and the molecules with the higher ex¬ 
citation energies are smoothly distributed over a similar range in column densities to their 
integrated intensities, their velocity dispersion distributions are typically concentrated at low 
values of the velocity dispersion. With the exception of HNCO which is far more abundant 
than these other molecules, we find that the emission predominantly arises from very narrow 
velocity features (a < 2 km s _1 ). Such narrow velocity dispersions have also been observed 
toward GO.253+0.016 on small spatial scales from other molecular tracers (i.e. N 2 H + (3-2); 
Kauffmann et al. 2013). The tail in these distributions to considerably higher velocity dis¬ 


persions (dotted contours), indicates positions within the clump where the moment analysis 
has overestimated the velocity dispersion due to multiple velocity components within the 
spectrum. This occurs in fewer than 1% of the spectra. 


The narrow range in dust column density over which the higher excitation energy molec¬ 
ular transitions are detected can be understood in terms of their critical densities. These 
molecules are typically detected when the gas has either a high temperature (T ffas >100 K) or 
a high density (n > n crit ), or both. Because these conditions are needed for both their forma¬ 
tion and excitation, the presence of these molecules can be used to pinpoint gas with these 
physical conditions. These molecular transitions all have critical densities of 10 1 2 * * 5 * — 10 7 cm 3 
(see Table [l]). Thus, while the gas temperature may be low toward the clump’s centre, the 
volume density may be very high. If the volume densities are sufficiently higher than their 
critical densities, the molecules may be excited regardless of the gas temperature. 

We see a clear trend that molecules with higher excitation energies and critical densities 


1 With the exception of HNCO, these plots were made using the ALMA only data and, thus, will not 
include the lower column densities traced on the large scale by the dust emission. Given the conditions 
needed to excite these transitions, large-scale emission is not expected. Indeed, their data do not show 
obvious effects that would indicate they suffer from the lack of the zero spacing information in a similar 
way to the data from the molecular transitions with lower excitation energies and the dust continuum very 
clearly do (i.e. large negative ‘bowls’ around bright emission). 

2 The distribution in measured velocities are narrower for HCN and HCO + compared to many of the 

other transitions because they include the large-scale emission from the single-dish data which, due to the 

‘intensity weighting’ of the moment calculation and in the case of widespread emission with a broad line 

profile, will tend to result in a similar calculated value for the intensity weighted velocity despite the emission 

being detected over a broad range. See Figure [9] for a comparison of these plots for both HCN and HCO + 

including and excluding the single dish data to see this effect. 
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peak toward the centre of the clump (compared to those with lower excitation energies and 
critical densities). The position-velocity diagrams also reveal that the molecular transitions 
with higher excitation energies and critical densities peak toward the centre of the cloud 
in both position and velocity (declination offset of ~ 0 , and velocity of ~ 20 km s -1 ; e.g. 
CH 3 CHO and NH 2 CN in Fig. g. These trends are exactly what is expected for a clump 
with a denser interior. 


5.2. Which molecules can be used to assess the dynamical state of the cores? 

While the molecular transitions detected by these observations all probe dense gas, 
they span a range in their critical densities and excitation energies, trace distinct chemical 
conditions, and often suffer from high optical depths. As such, they may often trace different 
material within GO.253+0.016 and may not necessarily correlate well with the optically thin 
dust continuum emission. 

One of our key goals is to identify star-forming dust cores within the clump, to measure 
their sizes, masses, and virial state. Determining these properties will allow us to ascertain 
the potential of GO.253+0.016 to form a cluster and the distribution of masses for its cores. 
In order to identify bound star-forming cores, we need to reliably associate molecular line 
emission with each dust continuum core in order derive their velocity dispersions. Ideally, 
we need to identify a molecular transition that is optically thin and is bright enough to be 
detected in most cores. 

Given the limited velocity resolution of our data, a detailed analysis of the cores’ virial 
states are challenging with the current dataj^] However, we can gauge, overall, which molec¬ 
ular transitions correlate with the dust emission. We find that emission from the molecular 
transitions with higher excitation energies typically provide a better morphological match to 
the dust column density compared to the transitions with lower excitation energies, the iso- 
topologues, or the shock tracers (compare the morphology of the dust emission to NH 2 CHO, 
CH 3 CHO, HNCO and H 2 CS in Fig. [3] and also their 2 -D correlation coefficients shown in 
Figure [ 4 ]). Moreover, since their velocity dispersions are typically narrow (~ 2 km s -1 ), these 
molecules may in fact be tracing individual high column density cores within the clump with 
very narrow velocity dispersions. 


3 The identification of discrete cores within GO.253+0.016 from the 3 mm dust continuum emission, an 
analysis of their virial state, and derivation of the core mass function will be presented in a future paper. 
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5.3. Is the gas structure set by turbulence? 


Supersonic turbulence is an important mechanism by which overdensities within molec¬ 
ular clouds may be generated (Vazquez-Semadeni 1994; Padoan et al. 1997). As such, un¬ 
derstanding the degree to which turbulence sets the initial structure within molecular clouds 
is an important part of understanding core and, ultimately, star formation. Recent work 
using the ALMA dust continuum emission presented here shows that the log-normal shape, 
dispersion, and mean of its column density probability distribution function (N-PDF) very 
closely matches the predictions of theoretical models of supersonic turbulence in gas given 
the high-pressure environment in which GO.253+0.016 is immersed (Rathborne et al. 2014b). 

These results also show that while, overall, the N- PDF shows a log-normal distribution, 
there is a small deviation from log-normal at the highest column densities, which indicates 
self-gravitating gas. This small deviation is coincident with the location of a known water 
maser (Lis & Carlstrom 1994) and, thus, pinpoints where star formation is beginning in the 
cloud. Assuming a dust temperature of 20 K, this small core (radius of 0.04 pc) has a mass 
of ~72M 0 and virial ratio, a v ; r , of < 1 (this is an upper limit since the associated molecular 
line emission is unresolved in velocity, a < 1.4km s -1 ; see Rathborne et al. 2014b for details). 
Since a v i r is close to unity, this core is gravitationally bound and unstable to collapse. 


The fact that the N-PDF for GO.253+0.016 matches well predictions from turbulent 
theory suggests that our current understanding of molecular cloud structure derived from 
the solar neighbourhood also holds in high-pressure environments (Rathborne et al. 2014b). 
This has important implications for understanding the initial gas structure for star formation 
in a range of different environments across the Universe (Kruijssen & Longmore 2013). 


To further determine the extent to which turbulence shapes the observed gas structure 
within GO.253+0.016, we have calculated the fractal dimension and the turbulent power spec¬ 
trum. These statistical methods are commonly used to describe the structure of molecular 
clouds in both simulations and observations. 


5.3.1. Fractal structure on small spatial scales 

An important characteristic of a protocluster clump, and one that is easily revealed 
by these data, is the distribution of its material on small spatial scales. If GO.253+0.016 
represents a clump that will give rise to a high-mass cluster, then the distribution of the 
material within it before star formation has commenced can provide important constraints 
for theories that describe cluster formation. As discussed previously, there have been a wide 
range of theoretical and numerical studies showing that the formation of high-mass stellar 
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clusters should proceed hierarchically, but thus far the degree of hierarchal structure within 
high-mass protoclusters has not been estimated. 


The fractal dimension (D f) can be used to characterise the structure within a molecular 
cloud. It is a simple tool that determines the degree to which the emission is smooth and 
centrally condensed or structured and hierarchical. Since observations of dust or molecular 
line emission reveal only the 2 dimensional projection of a 3 dimensional structure, care 
needs to be taken to accurately convert the measured projected fractal dimension (D p ) to 
the real three-dimensional fractal dimension (D f). This conversion is often assumed to be 


D/ = D p + 1 (Beech 1992); however, using simulations of cloud structure, Sanchez et al. 


(2005) derive an empirical relation to convert between D/ and D p (see their figure 8). Rather 
than a single conversion factor of unity for all values of D p , their simulations suggest that 
Dj may actually increase linearly as D p decreases. Previous observations derive projected 
fractal dimensions D p ~ 1.35 (see 
quoted D f for molecular clouds of 


Federrath et al. 2 

1009 

and references therein) leading to a 

~ 2.3 ( 

Elmegreei 

r 1997). 


To calculate the projected fractal dimension (D p ) for GO.253+0.016, we employ the 


perimeter-area method ( 

Mandelbrot] 1977 

which is often used to measure the fractal dimen- 

sion of interstellar gas clouds (e.g., 

Falgarone et al.||1991; 

Federrath et al. 2009; 

Sanchez et al. 


2005). Structures were defined as connected regions within a contour at a given intensity 


level. The area (A) was calculated as the sum of the number of pixels within the contour, 
the perimeter (P) was calculated as the number of pixels externally adjacent to the contour 
level. By varying the contour levels from the image minimum to its maximum, values for 
P and A were determined on many spatial scales. With a power-law fit of the form P oc A 
Dp / 2 , we can determine the perimeter-area fractal dimension D p . A smooth structure will 


have Dp = 1, while a very convoluted structure will have D„ = 2 (Federrath et al. 2009). 


Figure [12] shows the perimeter-area plots measured from contours of the dust and line 
integrated intensity emission. For this analysis we use, where possible, the data that includes 
the large-scale emission. Since this analysis is based on contours within an image, the effect 
of including the large-scale emission simply increases the number of contours to lower levels 
of the intensity which translates to an increase in the number of points plotted at large 
perimeters and areas. Their inclusion has little effect on the derived slope. 


In all panels the solid lines show two extreme values for D p (1 for smooth, centrally 
condensed emission, 2 for a highly fractal cloud). The dotted lines show linear fits to the 
data points: the derived values for D p range from 1.4-1.7, depending on the molecular tracer. 
For the dust continuum that includes the large-scale emission, we find a D p of 1.6 (plotted 
in the top left panel, marked as ALMA+SD). For comparison we also calculate the fractal 
dimension derived from only the large-scale dust continuum emission (in the same panel 
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marked as SD), and find D p to be ^1.2, consistent with its smooth distribution. 


Since the fractal dimension is best indicated by emission that is optically thin, i.e., from 
the dust emission and the molecular transitions with higher excitation energies, we select 
the derived D p from these molecules to characterise the emission from within GO.253+0.016 
(i.e., D p ~ 1.6). Using the empirical conversion from the two-dimensional fractal dimension 


Dp to the three-dimensional fractal dimension D/ from Sanchez et al. (2005), we find that 


D 


/ 


U 


for GO.253+0.016. Simulations of compressively driven turbulence in the ISM 
(Federrath et ah 2009) with Mach numbers, A4, of ~ 5.5, reveal typical values for D/ of 


2.3-2.6 (which correspond to D p ~ 1.3-1.5 using the empirical relation derived by Sanchez 
et al. 2005) . Since M. is ~ 30 for GO.253+0.016, it has a much higher level of turbulence 


compared to most other molecular clouds. Given this higher level of turbulence, it may 
be expected that the sonic length may be lower and the degree of fragmentation might be 
higher in GO.253+0.016 (i.e. a higher D p , and thus, lower D/ for higher Ai). While detailed 
simulations are necessary to determine the effect of A4 on the fractal dimension, our results 
are consistent with a marginally higher value for D p for GO.253+0.016 when compared to 
typical molecular clouds. This analysis suggests that the clump is highly fragmented on 
small spatial scales and its structure consistent with that produced by compressively driven 
turbulence. 


5.3.2. Turbulent power spectrum 

The power spectrum provides a measure of how turbulence is propagated within a 
molecular cloud and it identifies the preferential spatial scale, if any, at which the structure 
changes. On large scales a break in the power spectrum indicates the spatial scale at which 
turbulent energy is injected, on the small scales it can indicate the spatial scales at which 
either turbulence is dissipated or injected . Measuring the small-scale break point of the 
power spectrum can therefore reveal the spatial scale at which gravity overcomes the tur¬ 
bulent pressure or the spatial scale at which energy may be injected (i.e. from protostellar 
outflows) . 

We use the dust continuum image and several of the integrated intensity images to 
compute the power spectra. Because the turbulent power spectrum is a statistical measure 
of the energy on various spatial scales, it is important that the emission is well detected 
and well sampled across the clump. As such, we perform this analysis using HCN, HCO + , 
HNCO, and SiO because their data include the zero-spacing information and the emission 
is detected with high signal-to-noise ratios at all spatial scales across GO.253+0.016. 
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Figure [l3| shows the annular averages of each of the power spectra, P(k), as a function of 
the spatial frequency, k. The power spectra were generated using a Fast Fourier Transform of 
the images. In all panels the dotted vertical line marks the spatial frequency corresponding 
to the angular resolution of the data: the power spectra at spatial frequencies greater than 
this are not well sampled by these data. 


In cases where the power spectrum is close to a power law (i.e., for HCN, HCO + , HNCO, 
and SiO emission) we fit it with a single power law to derive a where, P(k) oc k~ a . In these 
cases, we find that a ~ 3. A power-law slope of 3 is predicted for a line intensity power 
spectrum of an optically thick medium (Lazarian & Pogosyan 2004; Burkhart et al. 2013). 
Thus, the derived power spectrum slope of ~ 3 for these molecular transitions are consistent 
with them being optically thick toward GO.253+0.016. 


In contrast, for the dust continuum emission which is optically thin, the power spectrum 
shows a clear break. While this break occurs at large spatial frequencies (i.e. small spatial 
scales), it is lies above the spatial frequencies corresponding to the angular resolution of these 
data (marked with the vertical dotted line). For this power spectrum we fit two distinct power 
laws and derive their indicies and the break point. For small values of k (i.e., large spatial 
scales) the derived power-law index is a ~ 2.9 compared to a ~ 7.5 for large values of k (i.e., 
small spatial scales). Values for a of ~ 2.8 are observed toward other star-forming regions 


such as Perseus, Taurus and Rosetta (Padoan et al. 2004). 


The trend of a steeper spectrum on smaller spatial scales implies that the higher column 
density material is concentrated on smaller spatial scales, while the lower column density 
material is distributed on larger spatial scales. This is consistent with the distribution of 
the emission in the dust continuum image. Since the effect of gravity on small scales will 
cause a steepening of the power-law slope at large spatial frequencies , the spatial scale at 
which gravity may dominate can be determined from the break point in the power spectrum. 
From the dust continuum emission power spectrum, we find that the break occurs at spatial 
frequencies corresponding to size scales of ~0.09pc. This size scale is comparable to the size 
of the small region of self-gravitating gas that was identified as the deviation at high column 


densities via the N-PDF (Rathborne et al. 2014b). As such, this lends support to the idea 


that gravity has overcome the internal pressure on these small spatial scales. 


For cores within molecular clouds to collapse into stars, they need to be sufficiently 
dense that their self-gravity overcomes the thermal pressure . The Jeans length, defined 
as the minimum length scale for gravitational fragmentation to occur, was determined for 
GO.253+0.016 using the relation 
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A j = OApc 


n{H 2 ) 


0.2 Arms -1 / V 10 3 cm -3 


- 0.5 


( 1 ) 


where c s is the sound speed and n(H 2 ) the volume density. While the measured dust temper- 


ature toward GO.253+0.016 is ~ 20K (Longmore et al. 

2012 

), recent observations indicate 

that the gas temperature is considerably higher, ~ 65 K ( 

Ao et al. 

2013 

), leading to the 


speculation that the gas may be heated by the dissipation of turbulent energy and/or cosmic 
rays rather than by photon heating. These observations are supported by SPH modelling of 
the dust and gas temperature distribution in GO.253+0.016 which suggest that the gas and 
dust temperatures are not coupled and that the high interstellar radiation held (ISRF) and 
cosmic ray ionisation rate (CRIR) in the CMZ are responsible for the discrepant dust and 
gas temperatures (Clark et al. 2013). 


Assuming a global gas temperature for GO.253+0.016 of 65K (Ao et al. 2013) the 
corresponding gas sound speed, c s , is 0.5 km s -1 . Using an average volume density for this 
material of ~ 10 5 cm -3 , we derive A j to be 0.10 pc. This value is comparable to the size scale 


derived from the break point in the power spectrum of the dust continuum emission and 
of the size of the self-gravitating core identified via the N-PDF. Extrapolating the global 
size-linewidth relation in Figure [l4] (top right panel) to smaller scales, suggests that the 
line-widths are subsonic on spatial scales of ~0.1pc (the sonic length is ~ 0.2pc). As 
such, since the sonic length exceeds the derived break point in the power spectrum, the 
thermal Jeans length is indeed the relevant quantity for describing the minimum size scale 
for gravitational fragmentation . An absence of a similar break in the power spectra derived 
from the molecular line emission is consistent with the transitions being optically thick and 
the dust continuum optically thin. Thus, the small, high-column density features detected 
in the optically thin dust continuum emission may have formed because gravity was able to 
overcome the thermal pressure on these small scales. 


6. Implications for the formation of young massive clusters 


Observations show that YMCs like the Arches are very centrally condensed, with stellar 
masses of >10 4 M 0 within their central OApc (Portegies Z wart et al.| 2010). If GO.253+0.016 
is the molecular cloud progenitor to a cluster like the Arches, then its internal structure can 
reveal the initial conditions for the formation of such high-mass clusters and whether or not 
the central concentration of stellar densities arises directly from the initial gas structure. 


To determine its internal structure, we have calculated the mass and volume density 
profiles of GO.253+0.016 as a function of radius (left panels of Fig.[l4|). We define each region 
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over which to calculate the mass and volume density using contours of column density; they 
were calculated on many spatial scales by varying the contour levels from the lowest to the 
highest values. An ‘effective radius’ was determined for each region by summing the number 
of pixels contained within that contour level and assuming that this area, A, forms a circle 
such that R e ff = \JA j". The mass within each region was then calculated by summing the 
column density, the volume density is calculated by dividing the mass by the volume of a 


sphere with a radius of R e ff (left panels of Fig. 14). 


Using the detected emission from the various molecular species, we have also calcu¬ 
lated the average velocity dispersion and virial ratio as a function of radius (right panels 
of Fig. 14). The virial ratio is tt v ; r , defined as 0 : vir = ha 2 R/GM, where a is the measured 


velocity dispersion, R the radius, G the gravitational constant, and M the mass. For these 
calculations the ‘effective radius’ was determined from the total number of pixels within the 
integrated intensity image for each of the molecular transitions. Within these regions the 
total mass, M, was determined from the dust column density summed over the same region, 
while cr was determined by averaging the intensity weighted velocity dispersions. 

The ALMA observations of the gas structure within GO.253+0.016 show that the mate¬ 
rial is distributed in a complex, filamentary network of structures that are highly fragmented. 


Its volume density profile (lower left panel of Fig. 14) shows that the average volume density 
increases with smaller radii, suggesting that the material is more concentrated toward the 
clump’s centre. We find that the volume density profile follows very well a power law for radii 
< 2 pc with an index of ~ —1.2. In contrast to what is expected for a Bonner-Ebert sphere 
and Plummer models that are often used to describe the volume density profiles within 
starless cores and globular clusters, the volume density profile for GO.253+0.016 shows no 
flattening at small radii. 

While the total mass of GO.253+0.016 is ~ 10 5 M Q , its mass profile shows that the 
majority of the mass originates from the low column density material in its outer edges, at 


large radii (upper left panel of Fig. 14). We find that the average velocity dispersion decreases 


on small spatial scales (upper panel of Fig. 14), similar to the results of Kauffmann et al. 


(2013). Moreover, we find that the derived virial ratios are > 1 for the molecular emission 


that is widespread (i.e. have large effective radii) but < 1 for the molecular species that are 


detected only over smaller regions (lower panel of Fig. 14). This trend of decreasing virial 
ratios as the radius decreases, implies that GO.253+0.016 is unbound at large radii, but 
bound at smaller radii. Thus, its outer envelope of lower density material may be unbound 
and expanding, while the inner denser region is bound and, in the absence of other supporting 
mechanisms, may be collapsing. 


The notion that the GO.253+0.016 may have a dense and collapsing interior with a low- 
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density surface that is expanding was posed recently by Bally et al. (2014) when discussing 
the nature and origin of the HCO + absorption filaments seen toward GO.253+0.016. This 
idea is supported by the observations of the large scale gas motions which show systemati¬ 
cally red-shifted optically thick emission compared to the optically thin and hot gas tracers 
(Rathborne et al. 2014a). In this scenario, the exterior layers of the cloud have a lower 
excitation temperature compared to the inner layers of the cloud, because the emission is 
sub-thermally excited (n<<n cr j t ). The expanding outer regions, coupled with a collapsing 
inner region, is consistent with the scenario in which GO.253+0.016 was formed recently as 
a result of a pericentre passage close to Sgr A* (Longmore et al.||2013; Kruijssen et al.||2014a) 


While a gravitationally bound and collapsing clump is consistent with what is expected 
in the very early stage of mass assembly within a protocluster, we find that there is no large 
concentration of mass in the centre of GO.253+0.016. Indeed, its total mass within a radius 
of 0.5 pc is only ~ 6 x 10 3 M Q : this is insufficient to form an Arches-like cluster by more 
than 2 orders of magnitude. Indeed, recent work characterising the structure within other 
high-mass, dense clumps in the CMZ show that the lack of a strongly peaked central mass 
distribution may be a general trend within these YMC-progenitor clouds ([Walker et al. sub¬ 


structure of the resulting cluster must evolve from this initial dispersed, hierarchical struc- 


ture to a smooth, centrally condensed configuration via a dynamical process (e.g. 

Allison 

et al. 2010; 

Moeckel & Bate 2010; Kruijssen et al. 2012; Parker & Dale 2013| Parker et al. 

2014 

. As such, any model for cluster formation that posits that stars form directly from an 


initially smooth, central concentration of gas with a mass similar to that of the final stellar 


cluster (e.g. 

Lada et al. 

1984; 

Geyer & Burkert 

2001; 

Boily & Kroupa 

2003; 

Goodwin & 

Bastian 2006; 

Baumgardt & Kroupa 

2007 

, is ruled out by these observations. 


mitted). Thus, if GO.253+0.016 is destined to form a high-mass (>10 4 M Q ) cluster, then the 


7. Conclusions 

Using observations from the new ALMA telescope we have conducted a statistical 
analysis of the small-scale structure and velocity within the high-mass protocluster clump 
GO.253+0.016. Its 3 mm dust continuum and molecular line emission reveal a complex net¬ 
work of emission features that have a complicated velocity structure. Filamentary structures, 
that are coherent in velocity, are seen in both emission and absorption. 

Widespread emission is detected from 17 different molecular transitions that probe a 
range in distinct physical and chemical conditions. Comparing their intensity, velocity dis¬ 
persion, and velocity fields, we find that the molecular ‘families’ trace the same type of 
material with the clump: transitions that are optically thick have broad velocity dispersions 
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and trace material on the large scales, while transitions that are optically thin have nar¬ 
row velocity dispersions and trace the densest material in the clump’s centre. Indeed, the 
dust column density is best traced by molecules with higher excitation energies and critical 
densities, consistent with a clump that has a denser interior. 

The emission is found to be highly fragmented (D/ ~ 1.8), consistent with structure 
that is produced by compressively driven turbulence. We find a clear break in the turbulent 
power spectra derived from the optically thin dust continuum emission at a spatial scale of 
~ 0.1 pc. This size scale is comparable to the Jeans length corresponding to the observed 
gas density and temperature, and also to the size of the small core of self-gravitating gas 
that was identified from the N-PDF. Thus, these structures might correspond to the spatial 
scale at which gravity has overcome the turbulent pressure. 

We speculate that GO.253+0.016 is on the verge of forming a cluster from hierarchical, 
filamentary structures that arise from a highly turbulent medium. Because the initial gas 
structure is fractal and hierarchical, its structure is quite different from the smooth, com¬ 
pact, very centrally condensed structure observed in high-mass clusters. Consequently, these 
observations rule out models of cluster formation that begin with such gas distributions. 
The implication is that if stars within YMCs form from hierarchical gas structures such as 
these observed in GO.253+0.016, then their distribution must change to a more centrally 
condensed configuration over a dynamical time scale. Thus, any theoretical study that aims 
to accurately characterise and understand the early phases of cluster formation ought to 
begin with an initial gas structure that is hierarchical rather than one that is smooth and 
centrally concentrated. 
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Table h Summary of the molecular transitions detected toward GO.253+0.016“. 
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The molecules are grouped according to the molecular ‘families’ and are ordered by increasing excitation energies. 
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Fig. 1.— Combined infrared and millimetre image of GO.253+0.016 (24 jum emission from 
the Spitzer Space Telescope shown as cyan traces the heated dust from the Galaxy; 3 mm 
emission from ALMA shown as red traces the dust from the clump’s interior) overlaid with 
contours of the single-dish JCMT 450 pm emission (contour levels are 30 to 90% of the peak, 
in steps of 20%). The scale bar marks a size scale of lpc. This cloud is so cold and dense 
it is seen as an extinction feature against the bright IR emission from the Galaxy. With 
the significant improvement in angular resolution and because ALMA sees through to the 
clump’s interior, we are now able to pinpoint the dust and gas within it and characterise its 
internal structure. 
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Fig. 2.— Images toward GO.253+0.016 showing both the significant improvement in angular 
resolution provided by ALMA and the recovery of the emission on the large spatial scales 
provided by the inclusion of the zero-spacing information. From left to right: Single dish, 
ALMA-only, and combined ALMA and single dish images. The upper row compares images 
of the 3 mm dust continuum emission (in units of mJybeam -1 ), while the lower row compares 
images of the HNCO integrated intensity emission (in units of mJybeam -1 km s -1 ). In all 
images the beam size is shown in the lower left corner. These (and all subsequent) images 
are shown in equatorial coordinates: the (0,0) offset position in R.A. and Dec is 17:46:09.59, 
-28:42:34.2 J2000. 
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Fig. 4.— 2-D cross-correlation coefficients for each tracer with respect to all others (the 
lower half of the figure shows the correlation coefficients via a colour scale, the upper half of 
the figure shows the derived values). For each particular tracer (labelled along the diagonal) 
the molecular transition for which it has the highest correlation coefficient can be found by 
searching for the darkest square either in the row to its left or in the column below (the 
corresponding values are mirrored in the upper portion of the figure). Overall, the dust 
continuum emission is not highly correlated with any gas tracer, but is most correlated with 
the morphology of the emission from the molecules with high excitation energies (NH 2 CHO 
in particular). As expected, the tracers within each molecular ‘family’ are most correlated 
with each other. 
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Fig. 5.— Example spectra selected from 6 positions across GO.253+0.016 (see Figure [3] and 
Table [2] for more details). The molecular transitions are ordered, from top to bottom, in 
increasing excitation energies. The dotted lines delineate the range in velocity over which 
the moments were computed (— 10 km s -1 <Vlsr < 70km s -1 ). Note that C 2 H has a second 
hyperfine component which is evident at blue shifted velocities. Because the transitions of 
NH 2 CHO and HC 5 N are close in frequency, evident in their spectra is the other transition. In 
each case the other emission line was excluded from the spectrum before the moment analysis 
was performed. Each spectrum is labelled indicating whether the large-scale emission was 
included (ALMA+SD) or not (ALMA only). Labels also indicated to which ‘family’ the 
molecular transitions belong. 
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Fig. 6.— Moment maps for HCN. Upper left: integrated intensity map (units of 
mJybeam _1 km s _1 ). Upper right: peak intensity map (units of mJybeam -1 ). Lower left: 
intensity weighted velocity field (units of km s _1 ). Lower right: intensity weighted dispersion 
(units of km s _1 ). Overlaid on the latter three maps are contours of the peak intensity (levels 
are 30, 50, 70 and 90% of the peak). In all images the angular resolution is shown in the 
lower left corner. 
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RA offset (decrees: J2000) 


Fig. 7.— Channel maps for HCN. Each panel shows the emission (in units of mJybeam -1 ) 
at the displayed velocity (in units of km s _1 ). The angular resolution is shown in the lower 
left corner of each panel. 
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Integrated intensity HCO + 
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Fig. 15.— Moment maps for HCO + . Upper left: integrated intensity map (units of 
mJybeam _1 km s _1 ). Upper right: peak intensity map (units of mJybeam -1 ). Lower left: 
intensity weighted velocity field (units of km s _1 ). Lower right: intensity weighted dispersion 
(units of km s _1 ). Overlaid on the latter three maps are contours of the peak intensity (levels 
are 30, 50, 70 and 90% of the peak). In all images the angular resolution is shown in the 
lower left corner. 
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Fig. 16.— Moment maps for C 2 H. Upper left: integrated intensity map (units of 
mJybeam _1 km s _1 ). Upper right: peak intensity map (units of mJybeam -1 ). Lower left: 
intensity weighted velocity field (units of km s _1 ). Lower right: intensity weighted dispersion 
(units of km s _1 ). Overlaid on the latter three maps are contours of the peak intensity (levels 
are 30, 50, 70 and 90% of the peak). In all images the angular resolution is shown in the 
lower left corner. 
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Fig. 17.— Moment maps for HCO. Upper left: integrated intensity map (units of 
mJybeam _1 km s -1 ). Upper right: peak intensity map (units of mJybeam -1 ). Lower left: 
intensity weighted velocity field (units of km s -1 ). Lower right: intensity weighted dispersion 
(units of km s -1 ). Overlaid on the latter three maps are contours of the peak intensity (levels 
are 70% of the peak). In all images the angular resolution is shown in the lower left corner. 
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RA offset (degrees; J2000) 


RA offset (degrees; J2000) 


Fig. 18.— Moment maps for SiO. Upper left: integrated intensity map (units of 
mJybeam _1 km s _1 ). Upper right: peak intensity map (units of mJybeam -1 ). Lower left: 
intensity weighted velocity field (units of km s _1 ). Lower right: intensity weighted dispersion 
(units of km s _1 ). Overlaid on the latter three maps are contours of the peak intensity (levels 
are 30, 50, 70 and 90% of the peak). In all images the angular resolution is shown in the 
lower left corner. 
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Fig. 19.— Moment maps for SO. Upper left: integrated intensity map (units of 
mJybeam _1 km s _1 ). Upper right: peak intensity map (units of mJybeam -1 ). Lower left: 
intensity weighted velocity field (units of km s _1 ). Lower right: intensity weighted dispersion 
(units of km s _1 ). Overlaid on the latter three maps are contours of the peak intensity (levels 
are 30, 50, 70 and 90% of the peak). In all images the angular resolution is shown in the 
lower left corner. 
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Fig. 20.— Moment maps for H 13 CN. Upper left: integrated intensity map (units of 
mJybeam _1 km s _1 ). Upper right: peak intensity map (units of mJybeam -1 ). Lower left: 
intensity weighted velocity field (units of km s _1 ). Lower right: intensity weighted dispersion 
(units of km s _1 ). Overlaid on the latter three maps are contours of the peak intensity (levels 
are 30, 50, 70 and 90% of the peak). In all images the angular resolution is shown in the 
lower left corner. 
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Integrated intensity H13CQ+ 
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Fig. 21.— Moment maps for H 13 CO + . Upper left: integrated intensity map (units of 
mJybeam _1 km s _1 ). Upper right: peak intensity map (units of mJybeam -1 ). Lower left: 
intensity weighted velocity field (units of km s _1 ). Lower right: intensity weighted dispersion 
(units of km s _1 ). Overlaid on the latter three maps are contours of the peak intensity (levels 
are 50, 70 and 90% of the peak). In all images the angular resolution is shown in the lower 
left corner. 
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RA offset (degrees; J2000) 
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Fig. 22.— Moment maps for HN 13 C. Upper left: integrated intensity map (units of 
mJybeam _1 km s _1 ). Upper right: peak intensity map (units of mJybeam -1 ). Lower left: 
intensity weighted velocity field (units of km s _1 ). Lower right: intensity weighted dispersion 
(units of km s _1 ). Overlaid on the latter three maps are contours of the peak intensity (levels 
are 30, 50, 70 and 90% of the peak). In all images the angular resolution is shown in the 
lower left corner. 
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Integrated intensity HCC13CN 



0.02 0.01 0 - 0.01 - 0.02 


RA offset (degrees; J2000) 


o 

H 

Q 



o 

CO 


o 


o 

C\2 


o 


o 

C\2 

I 


RA offset (degrees; J2000) 


Peak intensity HCC13CN 


o 

H 

Q 



0.02 0.01 0 - 0.01 - 0.02 
RA offset (degrees; J2000) 

Velocity dispersion HCC13CN 


o 

H 

Q 



0.02 0.01 0 - 0.01 - 0.02 
RA offset (degrees; J2000) 


Fig. 23.— Moment maps for HCC 13 CN. Upper left: integrated intensity map (units of 
mJybeam _1 km s _1 ). Upper right: peak intensity map (units of mJybeam -1 ). Lower left: 
intensity weighted velocity field (units of km s _1 ). Lower right: intensity weighted dispersion 
(units of km s -1 ). Overlaid on the latter three maps are contours of the peak intensity (levels 
are 30 and 70% of the peak). In all images the angular resolution is shown in the lower left 


corner. 
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RA offset (degrees; J2000) 
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Fig. 24.— Moment maps for HNCO. Upper left: integrated intensity map (units of 
mJybeam _1 km s _1 ). Upper right: peak intensity map (units of mJybeam -1 ). Lower left: 
intensity weighted velocity field (units of km s _1 ). Lower right: intensity weighted dispersion 
(units of km s _1 ). Overlaid on the latter three maps are contours of the peak intensity (levels 
are 30, 50, 70 and 90% of the peak). In all images the angular resolution is shown in the 
lower left corner. 
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Fig. 25.— Moment maps for CH 3 SH. Upper left: integrated intensity map (units of 
mJybeam _1 km s _1 ). Upper right: peak intensity map (units of mJybeam -1 ). Lower left: 
intensity weighted velocity field (units of km s _1 ). Lower right: intensity weighted dispersion 
(units of km s _1 ). Overlaid on the latter three maps are contours of the peak intensity (levels 
are 30 and 70% of the peak). In all images the angular resolution is shown in the lower left 


corner. 
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Integrated intensity NH2CHQ 
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Fig. 26.— Moment maps for NH 2 CHO. Upper left: integrated intensity map (units of 
mJybeam _1 km s -1 ). Upper right: peak intensity map (units of mJybeam -1 ). Lower left: 
intensity weighted velocity field (units of km s -1 ). Lower right: intensity weighted dispersion 
(units of km s -1 ). Overlaid on the latter three maps are contours of the peak intensity (levels 
are 50, 70 and 90% of the peak). In all images the angular resolution is shown in the lower 
left corner. 
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Fig. 27.— Moment maps for CH 3 CH0. Upper left: integrated intensity map (units of 
mJybeam _1 km s _1 ). Upper right: peak intensity map (units of mJybeam -1 ). Lower left: 
intensity weighted velocity field (units of km s _1 ). Lower right: intensity weighted dispersion 
(units of km s _1 ). Overlaid on the latter three maps are contours of the peak intensity (levels 
are 30, 50, 70 and 90% of the peak). In all images the angular resolution is shown in the 
lower left corner. 
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Integrated intensity H2CS 



'I I I I I I 1 I I I I 1 I I I I I I I I H _ 

0.02 0.01 0 - 0.01 - 0.02 


RA offset (degrees; J2000) 

Velocity field H2CS 



0.02 0.01 0 - 0.01 - 0.02 


RA offset (degrees; J2000) 


Peak intensity H2CS 




o 

0.02 0.01 0 - 0.01 - 0.02 
RA offset (degrees; J2000) 


Velocity dispersion H2CS 



0.02 0.01 0 - 0.01 - 0.02 
RA offset (degrees; J2000) 


Fig. 28.— Moment maps for H 2 CS. Upper left: integrated intensity map (units of 
mJybeam _1 km s _1 ). Upper right: peak intensity map (units of mJybeam -1 ). Lower left: 
intensity weighted velocity field (units of km s _1 ). Lower right: intensity weighted dispersion 
(units of km s _1 ). Overlaid on the latter three maps are contours of the peak intensity (levels 
are 30, 50, 70 and 90% of the peak). In all images the angular resolution is shown in the 
lower left corner. 

















- 59 - 


Integrated intensity NH2CN 


Peak intensity NH2CN 


u 

p 

p 



0.02 0.01 0 - 0.01 - 0.02 
RA offset (degrees; J2000) 


Velocity field NH2CN 


o 

p 

Q 



0.02 0.01 0 - 0.01 - 0.02 
RA offset (degrees; J2000) 


o 

p 

Q 



0.02 0.01 0 - 0.01 - 0.02 
RA offset (degrees; J2000) 

Velocity dispersion NH2CN 


o 

p 

Q 





0.02 0.01 0 - 0.01 - 0.02 
RA offset (degrees; J2000) 


Fig. 29.— Moment maps for NH 2 CN. Upper left: integrated intensity map (units of 
mJybeam _1 km s _1 ). Upper right: peak intensity map (units of mJybeam -1 ). Lower left: 
intensity weighted velocity field (units of km s _1 ). Lower right: intensity weighted dispersion 
(units of km s _1 ). Overlaid on the latter three maps are contours of the peak intensity (levels 
are 30 and 70% of the peak). In all images the angular resolution is shown in the lower left 


corner. 
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Fig. 30.— Moment maps for HC 5 N. Upper left: integrated intensity map (units of 
mJybeam _1 km s _1 ). Upper right: peak intensity map (units of mJybeam -1 ). Lower left: 
intensity weighted velocity field (units of km s _1 ). Lower right: intensity weighted dispersion 
(units of km s _1 ). Overlaid on the latter three maps are contours of the peak intensity (levels 
are 50, 70 and 90% of the peak). In all images the angular resolution is shown in the lower 
left corner. 

















DEC offset (degrees; J2000) 
- 0.02 0 0.02 


61 


-9.6 

■r i i i i i i i i i i i i i i i i i i i + 

-6.2 

■r i i i i i i i i i i i i i i i i i i i + 

-2.8 

-r , , , , i , , , , i , , , , i , , , , + 

0.6 

-r i i i i i i i i i i i i i i i i i i i * 

4.0 

r i i i i i i i i i i i i i i i i i i 11 

7.4 

... 

10.8 

. .. 

14.2 

. .. 

17.6 

r . .. 

21.0 

J/Kl •* " 

. .. 

24.4 

; i \ 

% 

... 

27.8 

f ; 

•Z^i, 

s 

4 

r , r , , 1 , , , , 1 , , , , 1 , , , , f 

^-1- 1 - 1- 1 -1- 3 

C0 ; 

r ‘ . ,r .‘ bo i 

: i 

34.6 

: 

- i^K - 

. .. , i iT, .t- 

.wM 

- co 

- CO 

41.4 

44.8 

48.2 

H 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 t 

51.6 

55.0 

; A 

... 

,,, 11. ,,,, i 

JM ' 

i 

* ,,, , i... ,,,, + 

* 

r . .. 

^ i 1 i 1 

: " 

58.4 

61.8 

65.2 

H 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 t 

68.6 

72.0 


. .. . 1 i i i '.. 

Hk 

. .. i i . ..+ 

-Tiiii .. 

■r .. 


0.02 0.01 0 - 0 . 01 - 0.02 
RA offset (decrees: J2000) 


Fig. 31.— Channel maps for HCO + . Each panel shows the emission (in units of mJy beam -1 ) 
at the displayed velocity (in units of km s -1 ). The angular resolution is shown in the lower 
left corner of each panel. 
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Fig. 32.— Channel maps for C 2 H. Each panel shows the emission (in units of mJybeam -1 ) 
at the displayed velocity (in units of km s -1 ). The angular resolution is shown in the lower 
left corner of each panel. 
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Fig. 33.— Channel maps for HCO. Each panel shows the emission (in units of mJy beam -1 ) 
at the displayed velocity (in units of km s -1 ). The angular resolution is shown in the lower 
left corner of each panel. 
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Fig. 34.— Channel maps for SiO. Each panel shows the emission (in units of mJybeam -1 ) 
at the displayed velocity (in units of km s -1 ). The angular resolution is shown in the lower 
left corner of each panel. 
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Fig. 35.— Channel maps for SO. Each panel shows the emission (in units of mJybeam -1 ) 
at the displayed velocity (in units of km s -1 ). The angular resolution is shown in the lower 
left corner of each panel. 

























- 66 - 



. .. . . . ^ t .I'ji?*',! t; t; » jl - 


0.02 0.01 0 - 0 . 01 - 0.02 
RA offset (degrees: J2000) 

Fig. 36.— Channel maps for H 13 CN. Each panel shows the emission (in units of mJy beam -1 ) 
at the displayed velocity (in units of km s -1 ). The angular resolution is shown in the lower 
left corner of each panel. 
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Fig. 37.— Channel maps for H 13 CO + . Each panel shows the emission (in units of 
mJy beam -1 ) at the displayed velocity (in units of km s -1 ). The angular resolution is shown 
in the lower left corner of each panel. 
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RA offset (decrees: J2000) 


Fig. 38.— Channel maps for HN 13 C. Each panel shows the emission (in units of mJy beam -1 ) 
at the displayed velocity (in units of km s -1 ). The angular resolution is shown in the lower 
left corner of each panel. 
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Fig. 39.— Channel maps for HCC 13 CN. Each panel shows the emission (in units of 
mJy beam -1 ) at the displayed velocity (in units of km s -1 ). The angular resolution is shown 
in the lower left corner of each panel. 
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Fig. 40.— Channel maps for HNCO. Each panel shows the emission (in units of mJy beam -1 ) 
at the displayed velocity (in units of km s -1 ). The angular resolution is shown in the lower 
left corner of each panel. 
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RA offset (decrees: J2000) 


Fig. 41.— Channel maps for CH 3 SH. Each panel shows the emission (in units of 
mJy beam -1 ) at the displayed velocity (in units of km s -1 ). The angular resolution is shown 
in the lower left corner of each panel. 

























72 



0.02 0.01 0 - 0 . 01 - 0.02 


o 

o 


C\2 

o 

o 


o 

o 


o 


RA offset (decrees: J2000) 


Fig. 42.— Channel maps for NH 2 CH0. Each panel shows the emission (in units of 
mJy beam -1 ) at the displayed velocity (in units of km s -1 ). The angular resolution is shown 
in the lower left corner of each panel. 
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Fig. 43.— Channel maps for CH 3 CH0. Each panel shows the emission (in units of 
mJy beam -1 ) at the displayed velocity (in units of km s -1 ). The angular resolution is shown 
in the lower left corner of each panel. 
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Fig. 44.— Channel maps for H 2 CS. Each panel shows the emission (in units of mJy beam -1 ) 
at the displayed velocity (in units of km s -1 ). The angular resolution is shown in the lower 
left corner of each panel. 
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Fig. 45.— Channel maps for NH 2 CN. Each panel shows the emission (in units of 
mJy beam -1 ) at the displayed velocity (in units of km s -1 ). The angular resolution is shown 
in the lower left corner of each panel. 
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Fig. 46.— Channel maps for HC 5 N. Each panel shows the emission (in units of mJy beam -1 ) 
at the displayed velocity (in units of km s -1 ). The angular resolution is shown in the lower 
left corner of each panel. 

























